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Abstract 

In 1982, Gessel showed that the Apery numbers associated to the ir¬ 
rationality of ("(S) satisfy Lucas congruences. Our main result is to prove 
corresponding congruences for all sporadic Apery-like sequences. In sev¬ 
eral cases, we are able to employ approaches due to McIntosh, Samol-van 
Straten and Rowland-Yassawi to establish these congruences. However, 
for the sequences often labeled Sis and {r}) we require a finer analysis. 

As an application, we investigate modulo which numbers these se¬ 
quences are periodic. In particular, we show that the Almkvist-Zudilin 
numbers are periodic modulo 8, a special property which they share with 
the Apery numbers. We also investigate primes which do not divide any 
term of a given Apery-like sequence. 


1 Introduction 


In his surprising proof [Ape79], [Poo79] of the irrationality of C(3), R. Apery 
introduced the sequence 


A{n) 



( 1 ) 


which has since been referred to as the Apery sequence. It was shown by I. Gessel 
[Ges82, Theorem 1] that, for any prime p, these numbers satisfy the Lucas 
congruences 

Ain) = A{nQ)A{ni) ■ ■ ■ Airij.) (modp), (2) 

where n = no -I- nip -I- • • • -|- is the expansion of n in base p. Initial work of 
F. Beukers [Beu02] and D. Zagier [Zag09], which was extended by G. Almkvist, 
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W. Zudilin [AZ06] and S. Cooper [Cool2], has complemented the Apery num¬ 
bers with a, conjecturally finite, set of sequences, known as Apery-like, which 
share (or are believed to share) many of the remarkable properties of the Apery 
numbers, such as connections to modular forms [SB85], [Beu87], [AOOO] or 
supercongruences [Beu85], [Cos88], [CCSIO], [OSll], [OS13], [OSS14]. After 
briefly reviewing Apery-like sequences in Section 2, we prove in Sections 3 and 
4 our main result that all of these sequences also satisfy the Lucas congruences 
(2). For all but two of the sequences, we establish these congruences in Section 3 
by extending a general approach provided by R. McIntosh [McI92]. The main 
difficulty, however, lies in establishing these congruences for the sequence ( 77 ). 
For this sequence, and to a lesser extent for the sequence sig, we require a much 
finer analysis, which is given separately in Section 4. 

In the approaches of Gessel and McIntosh, binomial sums, like (1), are used 
to derive Lucas congruences. Other known approaches to proving Lucas congru¬ 
ences for a sequence C{n) are based on expressing C{n) as the constant terms of 
powers of a Laurent polynomial or as the diagonal coefficients of a multivariate 
algebraic function. However, neither of these approaches is known to apply, 
for instance, to the sequence (ry). In the first approach, one seeks a Laurent 
polynomial A(a:) = A(a:i,..., Xd) such that C{n) is the constant term of A(a;). 
In that case, we write C{n) = ct A{x)'^ for brevity. If the Newton polyhedron 
of A(a;) has the origin as its only interior integral point, the results of K. Samol 
and D. van Straten [SvS09] (see also [MV13]) apply to show that C{n) satisfies 
the Dwork congruences 

C(jAm + n)C{\n/p\) = C(jA~^m + \ n/p\)C{n) (mod p’') (3) 


for all primes p and all integers m,n ^ 0, r ^ I. The case r = 1 of these 
congruences is equivalent to the Lucas congruences (2) for the sequence C{n). 
For instance, in the case of the Apery numbers (1), we have [StrI4, Remark 1.4] 


from which one may conclude that the Apery numbers satisfy the congruences 
(3), generalizing (2). Similarly, for the sequence (ry), one may derive from the 
binomial sum (22), using G. Egorychev’s method of coefficients [Ego84], that 
its nth term is given by ct A(a:, y, z)”, where 


A{x,y,z) 



1 


xy{l + z)® 


{I + x){l + y){l + z)'^ 


However, A(a:, y, z) is not a Laurent polynomial, and it is unclear if and how one 
could express the sequence (rf) as constant terms of powers of an appropriate 
Laurent polynomial. As a second general approach, E. Rowland and R. Yassawi 
[RY13] show that Lucas congruences hold for a certain class of sequences that 
can be represented as the diagonal Taylor coefficients of where s ^ 

1 is an integer and Q{x) € 'L[x\ is a multivariate polynomial. Again, while 
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such representations are known for some Apery-like sequences, see, for instance, 
[Strl4], no suitable representations are available for the sequences (rj) or sis- 
It was conjectured by S. Chowla, J. Cowles and M. Cowles [CCC80] and 
subsequently proven by I. Gessel [Ges82] that 


A{n) 


1, if n is even, 
5, if n is odd, 


(mod 8). 


(4) 


The congruences (4) show that the Apery numbers are periodic modulo 8, 
and it was recently demonstrated by E. Rowland and R. Yassawi [RY13] that 
they are not eventually periodic modulo 16, thus answering a question of Gessel. 
The Apery numbers are also periodic modulo 3 (see (47)) and their values 
modulo 9 are characterized by an extension of the Lucas congruences [Ges82]; 
see also the recent generalizations [KM15] of G. Krattenthaler and T. Miiller, 
who characterize generalized Apery numbers modulo 9. As an application of the 
Lucas congruences established in Sections 3 and 4, we address in Section 5 the 
natural question to which extent results like (4) are true for Apery-like numbers 
in general. In particular, we show in Theorem 5.3 that the Almkvist-Zudilin 
numbers are periodic modulo 8 as well. 

The primes 2, 3, 7,13, 23, 29,43,47,... do not divide any Apery number A{n), 
and E. Rowland and R. Yassawi [RY13] pose the question whether there are in¬ 
finitely many such primes. While this question remains open, we offer numerical 
and heuristic evidence that a positive proportion of the primes, namely, about 
^ 0.6065, do not divide any Apery number. In Section 6, we investigate 
the analogous question for other Apery-like numbers, and prove that Cooper’s 
sporadic sequences [Cool2] behave markedly differently. Indeed, for any given 
prime p, a fixed proportion of the last of the first p terms of these sequences 
is divisible by p. In the case of sums of powers of binomial coefficients, such a 
result has been proven by N. Calkin [Cal98]. 


2 Review of Apery-like numbers 

Along with the Apery numbers A(n), defined in (1), R. Apery also introduced 
the sequence 

which allowed him to (re)prove the irrationality of C(2). This sequence is the 
solution of the three-term recursion 

(n -I- l)^M„_|_i = {an^ + an + b)un — cn^w„_i, (5) 

with the choice of parameters (a, b, c) = (11, 3, —1) and initial conditions u_i = 
0, uq = 1. Because we divide by (n -|- 1)^ at each step, it is not to be expected 
that the recursion (5) should have an integer solution. Inspired by F. Beukers 
[Beu02], D. Zagier [Zag09] conducted a systematic search for other choices of 
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the parameters (a,b,c) for which the solution to (5), with initial conditions 
U-i = 0, uo = 1, is integral. After normalizing, and apart from degenerate 
cases, he discovered four hypergeometric, four Legendrian as well as six sporadic 
solutions. It is still open whether further solutions exist or even that there should 
be only finitely many solutions. The six sporadic solutions are reproduced in 
Table 1. Note that each binomial sum included in this table certifies that the 
corresponding sequence indeed consists of integers. 


(a, b, c) 

[Zag09] 

[AvSZll] 

A{n) 

(7,2,-8) 

A 

(a) 


(11,3,-1) 

D 

(b) 

?(:)’("r) 

(10,3,9) 

C 

(c) 

?©■(“) 

(12,4,32) 

E 

(d) 

/ 2k\ / 2(n — k)\ 

(9,3,27) 

B 

(f) 


(17,6,72) 

F 

(g) 

|;,-„-8.-(;)(‘)’ 


Table 1: The six sporadic solutions of (5) 


Similarly, the Apery numbers A{n), defined in (1), are the solution of the 
three-term recurrence 


(n -I- = (2n -I- l){an^ + an + b)un — n{cn^ + d)un-i, (6) 


with the choice of parameters (a, 6, c, d) = (17,5,1,0) and initial conditions 
U-i = 0, uq = 1. Systematic computer searches for further integer solutions 
have been performed by G. Almkvist and W. Zudilin [AZ06] in the case d = 
0 and, more recently, by S. Cooper [Cool2], who introduced the additional 
parameter d. As in the case of (5), apart from degenerate cases, only finitely 
many sequences have been discovered. In the case d = 0, there are again 
six sporadic sequences, which are recorded in Table 2. Moreover, by general 
principles (see [Cool2, Eq. (17)]), each of the sequences in Table 1 times 
is an integer solution of (6) with d yf 0. Apart from such expected solutions, 
Cooper also found three additional sporadic solutions, including 


sis{n) 



- 3fc - 1^ 


^2n — 


(7) 
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for n ^ 1, with sig(O) = 1, as well as s^ and sio, which are included in Table 2. 
Remarkably, these sequences are again connected to modular forms [Cool 2] (the 
subscript refers to the level) and satisfy supercongruences, which are proved in 
[OSS14]. Indeed, it was the corresponding modular forms and Ramanujan- 
type series for 1 /tt that led Cooper to study these sequences, and the binomial 
expressions for sy and sig were found subsequently by Zudilin (sequence sio was 
well-known before). 


(a, b, c, d) 

[AvSZll] 

A{n) 

(7,3,81,0) 

(S) 

B-)‘3"-"(i)("r)ii^ 

k \ / \ / 

(11,5,125,0) 

iv) 

|'-‘,‘(:)’((‘”-3r>(‘“3-;‘)) 

(10,4,64,0) 

(a) 

?(:)’(?)(t:(’) 

(12,4,16,0) 

(e) 


(9,3,-27,0) 

(C) 


(17,5,1,0) 

(7) 


(13,4,-27,3) 

S 7 

?(:)’(”r)(G) 

(6,2,-64,4) 

Sio 

?G)‘ 

(14,6,192,-12) 

S18 

defined in (7) 


Table 2: The sporadic solutions of (6) 


3 Lucas congruences 

It is a well-known and beautiful classical result of Lucas [Luc78] that the bino¬ 
mial coefficients satisfy the congruences 

C)^ (:)(::)■■■(:) 

where p is a prime and n^, respectively are the p-adic digits of n and k. That 
is, n = uq + riip -I- • • • -|- n^.p’’ and k = ko + kip -I- • • • -I- fc^p’' are the expansions 
of n and k in base p. Correspondingly, a sequence a{n) is said to satisfy Lucas 
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congruences, if the congruences 


a{n) = a{no)a{ni) ■ ■ ■ a{nr) (mod p) (9) 

hold for all primes p. It was shown by I. Gessel [Ges82, Theorem 1] that the 
Apery numbers A(n), defined in (1), satisfy Lucas congruences. E. Deutsch and 
B. Sagan [DS06, Theorem 5.9] show that the Lucas congruences (9) in fact hold 
for the family of generalized Apery sequences 



with r and s positive integers. This family includes the sequences (a), (b) from 
Table 1, and the sequences ( 7 ), sio from Table 2. The purpose of this section and 
Section 4 is to show that, in fact, all the Apery-like sequences in Tables 1 and 2 
satisfy the Lucas congruences (9). Using and extending the general framework 
provided by R. McIntosh [McI92, Theorem 6 ], which we review below, we are 
able to prove this claim for all of the sequences in the two tables, with the 
exception of the two sequences (rj) and Sig, for which we require a much finer 
analysis, which is given in Section 4. 

Theorem 3.1. Each of the sequences from Tables 1 and 2 satisfies the Lucas 
congruences (9). 

Remark 3.2. The Lucas congruences (9), in general, do not extend to prime 
powers. However, it is shown in [Ges82], and generalized in [KM15], that the 
Lucas congruences modulo 3 for the Apery numbers extend to hold modulo 9. 

On the other hand, numerical evidence suggests that all the Apery-like se¬ 
quences from Tables 1 and 2 in fact satisfy the Dwork congruences (3). While 
Theorem 3.1 proves the case r = 1 of these congruences, it would be desirable 
to establish the corresponding congruences modulo higher powers of primes. 

Following [McI92], we say that a function L : Z>q —)■ Z has the double Lucas 
property (DLP) if L{n, k) = 0, for k > n, and if 

L{n,k) = L(no,kQ)L{ni,ki) ■ ■ ■ L{nr,kr) (modp), (11) 


for every prime p. Here, as in ( 8 ), nt and ki are the p-adic digits of n and k, 
respectively. Equation ( 8 ) shows that the binomial coefficients (^) are a DLP 
function. More generally, it is shown in [McI92, Theorem 6 ] that, for positive 
integers tq, ri,..., r^. 


Lin, k) 




n -|- mk 
k 


( 12 ) 


is a DLP function. For instance, choosing the exponents as = 1, we find that 
the multinomial coefficient 

(n + mk)\ 
fc!"*+i(n- k)\ 


n + mk 

k,k,... ,k,n — k 
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is a DLP function for any integer m ^ 0. 

Suppose that L{n,k) is a DLP function and that G(n) and H{n) are LP 
functions, that is, the sequences G{n) and H{n) satisfy the Lucas congruences 
(9). Then, as shown in [McI92, Theorem 5], 


F{n) = Y. L{n, k)G{k)H{n — k) 


(13) 


fc =0 


is an LP function. Note that (12) and (13) combined are already sufficient 
to prove that the generalized Apery sequences, defined in (10), satisfy Lucas 
congruences. In order to apply this machinery more generally, and prove The¬ 
orem 3.1, our next results extend the repertoire of DLP functions. In fact, it 
turns out that we need a natural extension of the Lucas property to the case 
of three variables. We say that a function M : —>■ Z has the triple Lucas 

property (TLP) if M(n, k,j) = 0, for j > n, and if 


M(n,k,j) = M{no,ko,jo) ■ ■ ■ M{nr,kr,jr) (mod p), 


(14) 


for every prime p, where rzi, ki and ji are the p-adic digits of n, k and j, 
respectively. It is straightforward to prove the following analog of (13) for TLP 
functions. 

Lemma 3.3. If M{n,k,j) is a TLP function, then 

n 

L{n,k) = YFI{n,k,j) 
i=o 

satisfies the double Lucas congruences (11). In particular, if Lin, k) = 0, for 
k > n, then L{n, k) is a DLP function. 

Proof. Let p be a prime. It is enough to show that, for any nonnegative integers 
nQ,n',ko, k' such that uq < p and fco < P, 

L(nQ + n'p,ko + k'p) = L{no,ko)L{n',k') (modp). (15) 

Since the sum defining L{n, k) is naturally supported on j € {0,1,..., n}, we 
may extend it over all j G Z. Modulo p, we have 

L{n,k) = ^M(n,fc,j) 


p-i 

= YYl + fp) 

io=Oi'GZ 

ioGZ j'GZ 

= Lino,ko)L{n ,k'), 


which is what we had to prove. 


□ 
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Lemma 3.4. The function 


M{n,k,j) 



is a TLP function. 

Proof. Clearly, M{n,k,j) = 0, for j > n. In order to show that M{n,k,j) is a 
TLP function, we therefore need to show that, for any prime p. 


M{no+ n'p,ko + k'p,jo + fp) = M{no,ko,jo)M{n',k',j') (modp), (16) 


provided that 0 ^ ng, ko,jo < p and n', k',j' ^ 0. Observe that in the case jo > 
no both sides of the congruence (16) vanish because of the Lucas congruences (8) 
for the binomial coefficients. We may therefore proceed under the assumption 
that jo ^ no- 

Writing [x'^]f{x) for the coefficient of a;” in the polynomial f{x), we begin 
with the simple observation that 


= [x^]{l+xf+P 


Modulo p, we have 

(1 + x)'=+^ = (1 +x)'=«+^“(l +a;)('='+^')P = (1 + x)'=“+^«(l + xP)'='+^' (mod p). 

Since 0 ^ feg + jo < 2p, extracting the coefficient of x” = x”“(xP)” from this 
product results in the congruence 



Ao + jo\ fk' + j'\ ^ /fco + jo\ fk' + j'\ 
\ no )\ n' ) \no+pJ\n'-lJ 


(mod p). 


Note that, under our assumption that jo ^ ng, the second term on the right- 
hand side of this congruence vanishes (since uo + P ^ jo + P > jo + ko). This, 
along with (8), proves (16). □ 


Corollary 3.5. The function 

L{n, k) 



is a DLP function. 

Proof. Set j = fc in Lemma 3.4. 


Lemma 3.6. The function 




(3fc)! 

kP 


□ 


is a DLP function. 



Proof. Let p be a prime. As usual, we write n = ng+ri'p and k = kg + k'p where 
0 < uo < p and 0 ^ fco < p. In light of (8) and (13), the simple observation 



(17) 


demonstrates that the sequence of central binomial coefficients is an LP func¬ 
tion. We claim that 

(3A:)! _ fZk\ f2k\ 

~W “ 


is an LP function as well. From the Lucas congruences for the central binomials, 
that is 



C") 



(mod p). 


we observe that 
the congruences 


is divisible by p if 2fco ^ p. Hence, we only need to show 


(3fc)! ^ (3fco)! (3fc')! 


(mod p) 


(18) 


under the assumption that k^ <pl2. Note that 


= [x'^Ki+xr 

= [x'<^»{xPf]il+xf^°{l+xPf^' (modp) 

/3fco\ / 3 A:o \ / 3A:' \ / 3fco W 3fc' \ 

\ko)[k' J \ko+p) \k' -ij^ \ko + 2p) \k' - 2] ■ 


In the case /cq < p/2, we have /cq +P > 3fco, so that the last two terms on the 
right-hand side vanish. This proves (18). 

Next, we claim that 


/ n \ (3A;)! _ / no \ (S/cq)! \ (3^')! 

\3k) kP ~ \3ko) koP \3k') k'P 

By congruence (18), both sides vanish modulo p if 3fco 
if S/cq < Pi then the usual argument shows that 



(modp). (19) 

^ p. On the other hand. 



(mod p). 


In combination with (18), this proves (19). 

Finally, the congruences L{n,k) = L{no,ko)L{n',k'), that is 


^n—3k 



(3fc)! 

kP 



(3fco)! ^y_ 3 fc// n'\ 
koP \3P) 


(3fc0! 

k'P 


(mod p), (20) 


follow from Fermat’s little theorem and the fact that both sides vanish if 3fco > 
no or 3k' > n'. □ 
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We are now in a comfortable position to prove Theorem 3.1 for all but 
two of the sporadic Apery-like sequences. To show that sequences ( 77 ) and sig 
satisfy Lucas congruences as well requires considerable additional effort, and the 
corresponding proofs are given in Section 4. 

Proof of Theorem 3.1. Recall from (17) that the sequence of central binomial 
coefficients is an LP function. Further armed with (12) as well as Corollary 3.5 
and Lemma 3.6, the claimed Lucas congruences for the sequences (a), (b), (c), 
(d), (f), (a), (e), ( 7 ), Sio, sj follow from (13). It remains to consider the 
sequences (g), (d), (C) as well as ( 77 ) and Sis- 
Sequence (g) can be written as 


k=0 



where F{k) = ( 7 ) Franel numbers (sequence (a)), which we al¬ 

ready know to be an LP function. As a consequence of Fermat’s little theorem, 
the sequence aP is an LP function for any integer a. Hence, equation (13) 
applies to show that Ag(n) is an LP function. 

In order to see that sequence (i5) satisfies the Lucas congruences as well, it 
suffices to observe that L{n, k) = is almost a DLP function, that is, it 

satishes the congruences (11) but does not vanish for k > n. This is enough to 
conclude from Lemma 3.6 that 




n + k'\ (3fc)! 
k y “fc!^ 


is a DLP function. Since this is the summand of sequence (6), the desired Lucas 
congruences again follow from (13). 

On the other hand, for sequence (C), we observe that 


L{n,k) = 

3=0 


k + j 
n 


satishes the congruences (11) by Lemma 3.3 because the summand is a TLP 
function by Lemma 3.4. Hence, (^) L{n, k) is a DLP function. Writing se¬ 
quence (C) as 

” / \ 2 

AcH = E (r) 

the claimed congruences once more follow from (13). □ 


4 Proofs for the two remaining sequences 

The proof of the Lucas congruences in the previous section does not readily 
extend to the sequences ( 77 ) and sig from Table 2, because, in contrast to the 
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other cases, the known binomial sums for these sequences do not have the prop¬ 
erty that their summands satisfy the double Lucas property. Let us first note 
that the binomial sums for sis and sequence {ij), given in (7) and Table 2, can 
be simplified at the expense of working with binomial coefficients with negative 
entries. Namely, we have 


Sl8 


(n)=^(-l)‘ 


fc =0 


n\ / 2k\ f2{n — k)\ /2n — 3fc 


and 


A(n)=E(-l)' 


fc =0 


n — k 


An — 5k 
3n 


( 21 ) 


( 22 ) 


where, as usual, for any integer m ^ 0 and any number x, we define 


x{x — 1) • • • (x — TO -I- 1) 


For instance, the equivalence between (7) and (21) is a simple consequence of 
the fact that, for integers n ^ 0 and I = n — k, 


i-lf 


2n — 3k 
n 


= {-i) 


k-\-n 


—n + 3A: — 1 
n 


= (-i) 


I l^2n — 31 — 1 
n 


(23) 


For the first equality, we used that, for integers 6^0, 
a(a — 1) ■■■ (a — b + 1) 


bl 

^ , ^Nb (-a)(-Q + l)---(-Q + ^-1) ^ 
b\ 


—a + b — 1 
b 


(24) 


The following result generalizes the Lucas congruences for the sequence Sig(n). 

Theorem 4.1. Suppose that B(n, k) is a DLP function with the property that 
B{n, k) = B{n,n — k). Then, the sequence 


fc =0 \ ^ / 

is an LP function, that is, A{n) satisfy the Lucas congruences (9). 

Proof. Let p be a prime and let n ^ 0 be an integer. Write n = no + n'p and 
k = ko + k'p, where 0 Si no < p and 0 Si ko < p and n',k' are nonnegative 
integers. We have to show that 

A{n) = A{no)A{n') (modp). (25) 


In the sequel, we denote 

C{n, k) = {-lfB{n, k) ^ . 
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For fco ^ rro/3, we have 2no — S/cq ^ no ^ 0 and 2no — 3ko ^ 2no < ng + p. 
Hence, by the usual argument, we have 


^2n — Sk'j 


[x^°ixP f]{l + + xPf^’-^'^' (mod p) 



(mod p). 


Hence, we find that, when kg ^ ng/3, 

C(n,k) = C{no,kg)C{n\k') (modp). (26) 


For ng/3 < kg < 2ng/3, we have ng > 2ng — 3kg > 0. By the same argument as 
above, we find that 

f2n — 3k\ ^ 

f ^ ]=0 (modp), (27) 

and hence C{n^ k) = C{no, kg) = 0 modulo p. 

Finally, consider the case ng ^ 1 and 2no/3 ^ kg ^ ng. In that case, 
—p < —ng ^ 2ng — 3kg ^ 0 or, equivalently, 0 < 2no — 3A:o + p ^ p. Hence, we 

[^no( 2 .P)n'](i 2.)2no-3fco+P(^ ^p^2n'-3fc'-l 

/2no — 3kg +p\ /2n' — 3k' — 1 
V ng n' 

(2ng - 3kg\ /2n' - 3k' - 1 

\ ng J \ n' 

because, for any integers A, B and m such that 0 ^ m < p, 



have, modulo p, 


2n — 3k 
n 


/H + Bp\ 

\ m J 


^(H + Bp){A + Bp — 1) • • • (H + Bp — m + 1) 
ml 


—,A(A — 1) ■ ■ ■ (A — m-\-1) = 1 I (mod p). 
m\ \m-J 


(29) 


Set V = n' — k'. Applying (23) to the second binomial factor in (28), we find 
that 


^2n — 3fc^ 


(- 1 )”' 




(mod p). 


In combination with the assumed symmetry of B{n, fc), we therefore have that, 
when ng ^ 1 and 2ng/3 ^ kg ^ ng, 


C{n,k) = C{ng,kg)C{n',n'— k') (modp). (30) 
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We are now ready to combine all cases. First, suppose that ng ^ 1. Noting that 
k < n/3 implies k' < n'/3, and using (26), (27) and (30), we conclude that, 
modulo p, 


p—1 n' no n' 

fco—OA:'—0 kQ—Ok'—O 

[no/3j n' uq n' 

= Y 51 Y 

ko—O k'—O fco—r2no/3] fc'—0 

[no/3j n' no n 

= Yj C:{no,ko)'^C{n\k') + ^ C{no, kp) C{n', n' - k') 

fco—0 /c'—0 fco—r2no/3] k'—O 

[no/3j no n' 

Y, C(no,fco)+ 5Z 5]C(n',fc') 

fco—0 fco—|'2no/3] k'—O 

= A{no)A{n'), 


which is what we wanted to prove. The case rig = 0 is simpler, and we only 
have to use (26) to again conclude that (25) holds. □ 

Corollary 4.2. The sequence Si8(^T') satisfies the Lucas congruences (9). 


Proof. Recall from the discussion in Section 3 that 


B{n, k) 


/n\ /2(n — k)\ 

Uj vfcy V / 


is a DLP function. Obviously, B{n,k) = B{n,n — k). Hence, Theorem 4.1 
applies to show that Sis{n), in the form (21) satisfies the Lucas congruences 
(9). □ 

Next, we prove that the sequence ( 77 ), which corresponds to the choice a = 3 
in Theorem 4.3, satisfies Lucas congruences as well. 

Theorem 4.3. Let a G {1,3}. Then, the sequence 



(31) 


is an LP function, that is, A(n) satisfy the Lucas congruences (9). 

Proof. Let p be a prime and let n ^ 0 be an integer. As in the proof of 
Theorem 4.1, we write n = rig + n'p and k = kg + k'p, where 0 ^ ng < p and 
0 ^ kg < p and n', k' are nonnegative integers. Again, we have to show that 


A(n) = A(ng)A(n') (modp). (32) 
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Throughout the proof, let d = [3no/pJ. 

If fco ^ n-o/b, then 4no — 5fco ^ Suq ^ 0 and 4uo — 5^0 ^ 4no < Suq + p. 
Since d = [3uo/pJ j we thus have 0 ^ 3no — dp < p and 0 ^ 4no — S/cq — dp < 
(3no — dp) + p. Therefore, modulo p, 


An — 5k\ 
. 3n ) 




■x) 


4no —5fco —dp 


(1 


p\4n' —5A;'+d 


/4no — 5fco — dp\ /4n' — 5/c' + d\ 
\ 3no — dp / \ 3n' + d ) 

( 4no — 5ko'\ /4n' — 5k' + dX 
\ 3no ) \ 3n' + d / ’ 


where in the last step we used that, modulo p, 


Auq — 5ko — dp 
3nQ — dp 


Ano — 5ko — dp 
no - 5fco 


4no — bko 
no - 5ko 




which follows from (29) because 0 ^ no — 5ko < p. In particular, we have 


L^o/Sj n' 

E E(-i) 

ko—O k'—O 
[rio/Sj 

= E (-1)'“ 

fco=0 


4n — 5k 
3n 


no\^ ( 4no — 5A:c 


kf 


3nn 




k'—O 


/n'\ “ An' - 5k' + d\ 
\k') \ 3n' + d )’ 


(34) 


and we observe that, for d S {0,1}, 


A(n)=^(-I)^ 


fc=0 


4n — 5fc + d 
3n + d 


(35) 


To see this, note that the the sum of the k-th and (n — A:)-th term does not 
depend on the value of d e {0,1}. Indeed, using (24), Pascal’s relation and (24) 
again, we deduce that 


4n-5fc+l\ /4n-5(n-fc) + I 

3n + l /”*" V 3n + I 


4n — 5fc + 1 
3n + 1 
4n — 5fc + 1 
3n + 1 


4n — 5fc — 1 
3n + I 
4n — 5k 
3n + 1 


fAn- 

3n- 


4n — 5k\ f An — 5k — 1 

3n j V 3n 


An-5k\ /An-5{n- k) 

3n J ^ [ 3n 


5k 

- 1 


4n — 5fc — 1 
3n + 1 
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Next, suppose that no ^ 1 and 4no/5 ^ /cq ^ no- In that case, —p < —no ^ 
4no — 5ko ^ 0 or, equivalently, 0 < 4no — 5^o + P ^ P- Hence, we have, modulo 
P, 


4n — 5k\ 
. 3n ) 


+ a:)^"o-5fco+p^2^ _j_ ^p^4Ti'-5fc'-i 


4no — 5ko + p\ f ‘in' — 5k' — 1 


3no — dp 


3n' + d 


We rewrite the first binomial factor as follows, applying first (24) and then (29) 
twice, to find that, with Iq = ng — fcg, modulo p. 


4no — 5A:o + p 
3no — dp 


/_lNno+d [- 5Zo - (d + l)p - 1 

\ 3no — dp 


^ (- 1 ) 


no+d 


f 4no — 5Zo — dp — 1 
\ 3no — dp 


= (- 1 ) 


(_l)»o+d f 4no - 5Zo - dp - 1\ 

V no -5lo- 1 J 

(_l)»o+d/^4no-5Zo-l\ 

V no - 5^0 - 1 / 

/4no - 5Zo - 1^ 


ng+d I 


3no 


Here, we proceeded under the assumption that no — 5lo > 0. It is straightforward 
to check that the final congruence also holds when ng = 5lo, because then the 
binomial coefficients vanish modulo p. We conclude that, when ng ^ I and 
4no/5 ^ fcg ^ no, 


(-1)' 


4n — 5k 
3n 


^ (- 1 )' 


4no — 5^0 — I 
3no 


(- 1 ) 


k'+d 


4n' -5k' -1 
3n' + d 


In particular, we have 


(mod p). 


no n 


i: 

fco=|’4no/5] k'—O 


An — 5k 
3n 


no\“/4no - 5lo - 1 

lo J \ 3 no 




-\-d 


k'^0 


k' 


= Y. (-1)'“ 

fco=[4rao/5] 


“ /An' -5k' -1 


3n' + d 


k'—O 


“ ^4n' - 5A:' - I 
3n' + d 


(36) 


and we observe that, for integers d ^ 0, 


n / \ ^ 


k=0 


An — 5k — 1 
3n + d 




k=0 


An — 5k + d 
3n + d 
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because, by (24), 


(- 1 )^ 


An — hk + d 
in + d 




4n — 5(ri — fc) — 1 
3n + d 


Therefore, we can combine (34) and (36) into 


E 

fco—0 

or fco^4no/5 



^ A(no)E(-l)'' 

fc'=0 


/4n' - bk' + d\ 

\ 3n' + d J 


(mod p), 


(37) 


which holds for all 0 ^ ng < p (recall from the discussion at the beginning of 
this section that A{no), like sequence ( 77 ), can be represented as in Table 2). 

On the other hand, suppose that ng/S < ko < 4no/5. Set / = [(4ng — 
5A:g)/pJ. Since 0 < 4ng — b/cg < 3no < 3p, we have / G {0,1,2}. The usual 
arguments show that, modulo p, 


An — 5k\ 
. 3n ) 


|■2;3no-<ip^2;^’)^"'+‘^](l + x)'^'r>-o-5ko-fp _|_ ^p^4n'-5/c'+/ 


4ng — 5fcg — fp\ f An' — 5k' + f 


3rig — dp J \ 3n' 

/4ng — 5fcg\ / 4n' — 5k' + f 
\ 3no — dp ) \ 3n' + d 


(38) 


We are now in a position to begin piecing everything together. To do so, we 
consider individually the cases corresponding to the value of d S (0,1, 2}. 

First, suppose d = 0 or d = 1. Congruence (37) coupled with (35) implies 
that 


ko^no/5 or /co^4no/5 

To conclude the desired congruence (32), it therefore only remains to show that 


[4rio/5]-l n' 

E 

feo = l»io/5j+l k’—O 


An — 5k 
3n 


= 0 (mod p). 


(39) 


This is easily seen in the case d = 0, because then each term of this sum vanishes 
modulo p. Equivalently, for d = 0, (38) vanishes whenever no/5 < ko < 4ng/5 
(because 0 < 4ng — 5A;g — fp ^ 4ng — 5fcg < 3no). On the other hand, if d = 1, 
we claim that the sum (39) vanishes modulo p because the terms corresponding 
to {ko, k') and {ko, n' — k') cancel each other. To see that, observe first that, for 
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d = 1, (38) vanishes whenever no/5 < ko < 4no/5 and / = [(4no — 5ko)/p\ ^ 0 
(because 0 < 4no — 5ko — fp ^ 4no — 5ko—p < Zn^—p if / G {1,2}). Therefore, 
for the term corresponding to {kQ,k'), 


{-If 


An — 5k\ 
. 3n ) 


(-l)ko /^4no-5fco\ fc. 

^ ^ \3no-pJ^ ’ 


/ 4n' — 5fc'\ 
V 3n' + 1 ) 


(mod p), 


while, for the term corresponding to {ko, n' — k') with j = ko + {n' — k')p, 


(- 1 ) 


j ^ 4n - 5j 
3n 


= r_i')feo ( 4^^° - 5/co\ ( 4n' - 5(n' - k') 

~ ^ ’ \Zno-p ^ V Zn' + l 

V Zno-p 

- -(-!)'= (-dp), 


f_nfc'+i(^dn'-5fc'\ 

^ ^ \3n' + 1 J 


where we applied (24) for the second congruence. It is now immediate to see 
that the sum (39) indeed vanishes modulo p for d = 1. 

It remains to prove the Lucas congruences (32) in the case d = 2. Using 
(37), we have 


k'=0 \ / \ 


4n' - 5/fc' + 2 
3n' + 2 


M (mod p), 


where 


[4no/5]-l n' 

M--= E E(-i) 

fco=L"o/5j+l fc'=0 

Combining this congruence with the identity 


n\ [An — 5k 
3n 


k=0 ^ ^ 


/4n - 5A: + 2\ 
V 3n + 2 ) 


/An — 5k\ 

\3n + 2)\' 


which can be deduced along the same lines as (35), we find that 


A{n) = A{no)A{n') + A{no) E (“I)'"' 


k'=0 


“ /4ri' - 5k' 
3n' + 2 


+ M (modp). (40) 


We have, by (38), modulo p, 

r4rio/5]-l 


M= E (-1)''“ 

ko = [no/5\+l 

r4no/5l-l 

^ E (-1)'“ 

feo = L"o/5j+l 


no\ “ /4no - 5ko 

ko) V - 2p 

no\ “ /4no - 5ko 
ko) V - 2p 


n' 

i:(-u 


“ ^4n' — 5k' + / 
3n' + 2 


fe '=0 


“ /^n' - 5k'\ 
3n' + 2 /’ 
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where the last congruence is a consequence of the identity 



/4n — 5fc + 1\ f n\‘^ f An — 5k\ 

( 3„ + 2 U) (3.. + 2j 


(which follows from (24) and replacing k with n — k) and the fact that (38) 
vanishes for no/5 < ko < Ano/5 if / = 2. Using this value of M in (40), we find 
that the desired Lucas congruence (32) follows, if we can show that 


r4no/5]-l / \ a r-; \ 

AK)+ y: ^ »)s0 (modp)^ (41) 

fco=[no/5j+l \ c/ \ u V/ 

Note that, if ko ^ no/5, then, by (29) and (33), 

Mn, - 5t„\ ^ /4„„ - St. - 2p\ ^ /4n„ - 5t„\ 

\3no-2pJ V 3no-2p y V 3no y ^ ^ ' 

A similar argument, combined with (24), shows that the congruence (42) also 
holds if ko ^ 4no/5. We therefore find that (41) is equivalent to 



(mod p). 


The next lemma proves that this congruence indeed holds provided that a G 
{1,3}. □ 

Lemma 4.4. Let p be a prime, and a € {1, 2, 3}. Then we have, for all n such 
that 2p/3 ^ n < p, 



(mod p). 


Proof. To prove these congruences we employ N. Calkin’s technique [Cal98] for 
proving similar divisibility results for sums of powers of binomials (58). Denoting 
r = p — n, we have, by (24) and (29), 




p—r 


-E(-i) 

k—O 
p—r 


k (p — r\^ f Ap — Ar — 5k 


k—Q 

p—r 


k J \ p — 3r 

k — p + r — 1\"' f Ap — Ar — 5k 
p — 3r 


.^/fc + r-l\ /4p-4r- 5k\ 
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Clearly, 


k + r — l\ {k + l)(/c + 2) • • • (fc + r — 1) (fc + l)r-i 


(r- 1)! 


(r-1)! ’ 


(43) 


where {x)k = x{x -\- 1) ■■■ {x + k — 1) denotes the Pochhammer symbol (in 
particular, (a;)o = 1). Likewise, 

4p — 4r — 5fc\ (3p — r — 5A: + l)p- 3 r 
p—5r ) (p—3r)! 

Since (r — 1)! and {p — 3r)! are not divisible by p, we have to show that 


y^(fc + l)“_i(3p - r - 5k + l)p_ 3 r = 0 (mod p). 


(44) 


fc=0 


Since the polynomials {x)k, {x)k-i, ■ • •, {x)q form an integer basis for the space 
of all polynomials with integer coefficients and degree at most fc, there exist 
integers cq, Ci,..., cat with N = {a — l)(r — 1) + p — 3r so that 


N 


{k + l)“_i(3p - r - 5fe + 1)^-3^ = ^ Cj{k + r)j. 

j=o 


Then the left-hand side of (44) becomes 


N 


+ l)r-i'^Cj{k + r)j 
3=0 


fc=0 


= E 


N p—r 

'^Cj'^{k + l)r-l{k- 
j—0 k—0 

N p—r 

j—0 k—O 

^ (p - r + l)r+j 


3=0 


r + j 


where we used 

to evaluate 

p—r 


{x)k - {x- l)fe = k{x)k-i 


+ l)r+_;-l — E! 


{k l)r+j ik)r+j _ (P T Or+j 


r+j r+j 

The desired congruence therefore follows if we can show that 

(p - r + l)r+j 


r+j 


= 0 (mod p) 


(45) 


(46) 
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for all j = 0,1,... ,N. Since r > 0 and j ^ 0, the numerator {p — r + l)r+j is 
always divisible by p. The congruences (46) thus follow if r + j < p for all j, or, 
equivalently, r + N < p. Since 

r + = (a — l)(r — 1) + p — 2r, 

we have r + N < p if and only if 

(a — l)(r — 1) < 2r. 

Clearly, this inequality holds for all r ^ 1 if and only if a ^ 3. □ 

Remark 4.5. Numerical evidence suggests that the values a G {1,3} in Theo¬ 
rem 4.3 are the only choices for which the sequence (31) satisfies Lucas congru¬ 
ences. In light of Lemma 4.4, it is natural to ask if there are additional values 
of a and e, for which the sequence 



satisfies Lucas congruences. Empirically, this does not appear to be the case. In 
particular, for a = 2 this sequence does not satisfy Lucas congruences for either 
e = 0 or e = 1. 

5 Periodicity of residues 

The Apery numbers satisfy 

A(n) = (—1)” (mod 3), (47) 

and so are periodic modulo 3. As in the case of the congruences (4), which 
show that the Apery numbers are also periodic modulo 8, the congruences (47) 
were first conjectured in [CCC80] and then proven in [Ges82]. We say that 
a sequence C(n) is eventually periodic if there exists an integer M > 0 such 
that C{n + M) = C(n) for all sufficiently large n. An initial numerical search 
suggests that each sporadic Apery-like sequence listed in Tables 1 and 2 can 
only be eventually periodic modulo a prime p if p ^ 5. As an application of 
Theorem 3.1, we prove this claim next. 

Corollary 5.1. None of the sequenees from Tables 1 and 2 is eventually periodic 
modulo p for any prime p > 5. 

Proof. Gessel [Ges82] shows that, if a sequence C'(n) satisfies the Lucas con¬ 
gruences (9) modulo p and is eventually periodic modulo p, then C(n) = C(l)"’ 
modulo p for all n = 0,1,... ,p — 1. 

For instance, let C{n) be the Almkvist-Zudilin sequence (d). Then, C(l) = 
3, C{2) = 9 and C(3) = 3. Suppose C{n) was eventually periodic modulo p. 
Then p has to divide C(3) — C(l)^ = —24, which implies that p G (2, 3}. 
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In Table 3 we list, for each sequence, the primes dividing both C(2) — C(l)^ 
and (7(3) — (7(1)^. The fact, that all these primes are at most 5, proves our 
claim. □ 


(a) 

(b) 

(c) 

(d) 

(f) 

(g) 

(^) 

( 7 ) 

(«) 

(e) 

(C) 

( 7 ) 

(st) 

(sio) 

(sis) 

2,3 

2,5 

2,3 

2 

2,3 

2,3 

2,3 

2,5 

2,3 

2,3 

2,3 

2,3 

2 

2 

2,3 


Table 3: The primes dividing both (7(2) — (7(1)^ and (7(3) — (7(1)^, for each 
sequence C{n) from Tables 1 and 2. 


In fact, as another simple consequence of Theorem 3.1, we observe that the 
Apery-like sequences are in fact eventually periodic modulo each of the primes 
listed in Table 3. 


Corollary 5.2. Let C(n) be any sequence from Tables 1 and 2. 

• (7(n) = (7(1) (mod 2) for all n ^ 1. 

• (7(n) = (7(1) (mod 3) for all n ^ I if C{n) is one of (c), (/), {g), (S), 
(a), (e), (^), Sig, and C{n) = (—1)" (mod 3) for all 0 if C{n) is (a) 
or ( 7 ). 

• C(n) = 3" (mod 5) for all n ^ 0 if C{n) is {b), and C(n) = 0 (mod 5) 
for all n ^ 1 if C(n) is (ry). 

Proof. One can check that Table 3 does not change if we include only those 
primes p such that (7(n) — (7(1)"' is divisible by p for all n G {0,1, 2, 3,4}. For 
n = 0, this is trivial since (7(0) = 1. Therefore, in each of the cases considered 
here, we have 

(7(n) = (7(1)" (mod p) 

for all n G {0,1,... ,p — 1}. For any n ^ 0, let n = uq -I- nip -I- • • • -I- Urp’' be the 
p-adic expansion of n. Then, by Theorem 3.1, we have 

C(n) = C{no)C(ni) ■ ■ ■ C(ur) (mod p) 

= C'(l)"“+"i+-+"’- (mod p) 

= (7(1)" (modp). 


For the final congruence we used Fermat’s little theorem. All claimed congru¬ 
ences then follow from the specific initial values of (7(n) modulo p. □ 

More interestingly, the congruences (4) show that the Apery numbers (se¬ 
quence ( 7 )) are periodic modulo 8 . We offer the following corresponding result 
for the Almkvist-Zudilin sequence (S). 

Theorem 5.3. The Almkvist-Zudilin numbers 


Z{n) = ^(-1)'=3"-3'= 

k=0 



n + k\ (3fc)! 
n J kP 
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satisfy the congruences 


1 , if n is even, 
3, if n is odd, 


(mod 8 ). 


(48) 


Z{n) = 


Proof. It is shown in [Strl4] that the numbers (—l)"^(n) are the diagonal 
Taylor coefficients of the multivariate rational function 


F{xi,X2,X3,X4) 


1 

I — {xi + X2 + X3 + Xi) + 21X1X2X3X4 


(49) 


That is, if 


00 00 00 00 

F{xi,X2,X3,X4) = X! X! X! C(ni,n2,n3,n4)xf^X2'^x'tf^x1'^ 

ni=0 712—0 n3=0 714,-0 

is the Taylor expansion of the rational function F, then Z{n) = (—l)"C'(n, n, n, n). 

Given such a rational function as well as a reasonably small prime power p’’, 
Rowland and Yassawi [RY13] give an explicit algorithm for computing a finite 
state automaton, which produces the values of the diagonal coefficients modulo 
. In the present case, this finite state automaton for the values (—I)"Z(n) 
modulo 8 turns out to be the same automaton as the one for the Apery numbers 
modulo 8 . Hence, the congruences (48) follow from the congruences (4). We 
refer to [RY13] for details on finite state automata and the algorithm to construct 
them from a multivariate rational generating function. □ 


Empirically, Theorem 5.3 is the only other interesting set of congruences, 
apart from the congruences (4) , which demonstrates that an Apery-like sequence 
is periodic modulo a prime power. More precisely, numerical evidence suggests 
that none of the sequences in Tables 1 and 2 is eventually periodic modulo 
, for some r > 1, unless p = 2. Moreover, the only other instances modulo 
a power of 2 appear to be the following, less interesting, ones: sequences (d) 
and (a) are eventually periodic modulo 4 because all their terms, except the 
first, are divisible by 4; likewise, sequences (e) and S 7 are eventually periodic 
modulo 8 because all their terms, except the first, are divisible by 8 . We do 
not attempt to prove these claims here. We remark, however, that these claims 
can be established by the approach used in the proof of Theorem 5.3, provided 
that one is able to determine a computationally accessible analog of (49) for the 
sequence at hand. 


6 Primes not dividing Apery-like numbers 

Using the Lucas congruences proved in Theorem 3.1, it is straightforward to 
verify whether or not a given prime divides some Apery-like number. 

Example 6.1. The values of Apery numbers A(0), A(l),..., A( 6 ) modulo 7 are 
1,5, 3, 3, 3, 5,1. Since 7 does not divide A(0), A(l),..., A( 6 ), it follows from the 
Lucas congruences (9) that 7 does not divide any Apery number. o 
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Arguing as in Example 6.1, one finds that the primes 2, 3, 7,13,23,29,43,47,... 
do not divide any Apery number A{n). E. Rowland and R. Yassawi [RY13] pose 
the question whether there are infinitely many such primes. Table 4 records, 
for each sporadic Apery-like sequence, the primes below 100 which do not di¬ 
vide any of its terms, and the last column gives the proportion of primes below 
10^ with this property. Each Apery-like sequence is specified by its label from 
[AvSZll], which is also used in Tables 1 and 2. The alert reader will notice that 
Cooper’s sporadic sequences (the ones with d ^ 0 in Table 2) are missing from 
Table 4. That is because these sequences turn out to be divisible by all primes. 

A more precise result for these sequences is proved at the end of this section. 


(a) 

3, 11, 17, 19, 43, 83, 89, 97 

0.2994 

(b) 

2, 5, 13, 17, 29, 37, 41, 61, 73, 89 

0.2897 

(c) 

2, 7, 13, 37, 61, 73 

0.2962 

(d) 

3, 11, 17, 19, 43, 59, 73, 83, 89 

0.2815 

(f) 

2, 5, 13, 17, 29, 37, 41, 61, 73, 97 

0.2994 

(g) 

5, 11, 29, 31, 59, 79 

0.2929 

(^) 

2, 5, 7, 11, 13, 19, 29, 41, 47, 61, 67, 71, 73, 89, 97 

0.6192 

iv) 

2, 3, 17, 19, 23, 31, 47, 53, 61 

0.2897 

(a) 

3, 5, 13, 17, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 83, 89 

0.5989 

(e) 

3, 7, 13, 19, 23, 29, 31, 37, 43, 47, 61, 67, 73, 83, 89 

0.6037 

(C) 

2, 5, 7, 13, 17, 19, 29, 37, 43, 47, 59, 61, 67, 71, 83, 89 

0.6046 

( 7 ) 

2, 3, 7, 13, 23, 29, 43, 47, 53, 67, 71, 79, 83, 89 

0.6168 


Table 4: The primes below 100 not dividing Apery-like numbers (sequence indi¬ 
cated in first column using the labels from [AvSZll]) as well as the proportion 
of primes (in the last column) below 10,000 not dividing any term 

Example 6.1 shows that the first 7 values of the Apery numbers modulo 7 
are palindromic. Our next result, which was noticed by E. Rowland, shows that 
this is true for all primes. 

Lemma 6.2. For any •prime p, and integers n such that 0 ^ n < p, the Apery 
numbers A{n) satisfy the congruence 

A{n) = A{p — \ — n) (modp). (50) 

Proof. For n such that 0 ^ n < p, we employ (24) and (29) to arrive at 

^ fp - 1 - n\^ fp - I - n + k\‘^ 

= g( * j( . j 

^ gCt*") it) 


as claimed. 


□ 
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Theorem 3.1 and Lemma 6.2, considered together, suggest that « 

60.65% of the primes do not divide any Apery number. Indeed, let us make 
the empirical assumption that the values A(n) modulo p, for u = 0,1,..., (p — 
l)/ 2 , are independent and uniformly random. Since one of the values A{n) is 
congruent to 0 modulo p with probability 1 /p, it follows that the probability 
that p does not divide any of the (p + l )/2 first values is 

(l-p) ■ (51) 

By the Lucas congruences, shown in Theorem 3.1, and Lemma 6.2, p does not 
divide any of the (p+1)/2 first values if and only if p does not divide any Apery 
number. In the limit p—^oo, the proportion (51) becomes Observe that 

this empirical prediction matches the numerical data in Table 4 rather well. We 
therefore arrive at the following conjecture. 

Conjecture 6.3. The proportion of primes not dividing any Apery number 
A{n) is 



Figure 1: Proportion of primes (up to 10,000) not dividing the sequences (d), 
(a), (e), {(), ( 7 ), with the dotted line indicating The Apery sequence is 

plotted in blue. (We thank Arian Daneshvar for producing this plot.) 

While Lemma 6.2 does not hold for the other Apery-like numbers C{n) from 
Tables 1 and 2, we make the weaker observation that if a prime p > 5 divides 
C{n), where 0 ^ u < p, then p also divides C(p — 1 — u). We expect that this 
empirical observation can be proven in the spirit of the proof of Lemma 6.2, 
but do not pursue this theme further. We only note that it allows us to extend 
the heuristic leading to Conjecture 6.3 to the Apery-like sequences {6), (a), 
(e), (C) from Table 2. In other words, we conjecture that, for each of these 
sequences, the proportion of primes not dividing any of the terms is again 
Figure 1 visualizes some numerical evidence for this conjecture. On the other 
hand, for sequence ( 77 ) as well as the sequences from Table 1, the proportion of 
primes not dividing any of their terms appears to be about half of that, that is 
e-%V2 « 30.33%. 

To explain this extra factor of 1/2, we note that, for the Apery-like numbers 
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Stienstra and Beukers [SB85] proved that, modulo p, 


/p — 1 \ _ f 4a^ — 2p, if p = + b^, a odd, 

\^2 y~|o, ifp=3 (mod 4) 


(and conjectured that the congruence should hold modulo p^, which was later 
proved by Ahlgren and Ono [AOOO]). In particular, congruence (53) makes 
it explicit that every prime p = 3 (mod 4) divides a term of this Apery-like 
sequence. Note that, by a classical congruence of Gauss, the congruences (53) 
are equivalent, modulo p, to the congruences 


Ah 


(LiJ 


(b/2j 


0 , 


Vp/i\> 


if p = 1 (mod 4), 
otherwise, 


(54) 


which are valid for any prime p ^ 2. The more general result in [SB85] also 
includes the cases Aa and Ac- Similar divisibility results appear to hold for the 
other Apery-like numbers from Table 1, and it would be interesting to make 
these explicit. 

On the other hand, the extra factor of 1/2 in case of sequence (p) is explained 
by the following congruences, which resemble (54) remarkably well. 


Theorem 6.4. For any prime p 3, we have that, modulo p, 



Lp/3J V 

Lp/15J1 


if p = 1, 2,4, 8 (mod 15), 
otherwise. 


(55) 


Proof. Suppose that p = 2 (mod 3), and write p = 3n + 2. The congruence (55) 
can be checked directly for p = 2 and p = 5, and so we may assume p > 5 in 
the sequel. Applying (43) to the definition of sequence (p) in Table 2, we have 



4n — 5fc — 1\ 
3n ) 



- 5 fc) 3 „ (n 
(3n)! 


^4n — 5/c\ \ 

. 3n JJ 

5k + l) 3 n \ 
(3n)! J ■ 


(56) 


Since 3n = p — 2 and 0 ^ fc ^ n/5, the term 

jn - 5k)3n 
(3n)! 

is always divisible by p, unless n — 5k £ {1, 2} (because, otherwise, one of the 
p — 2 factors of (n — 5fc)3„ is divisible by p, while (3n)! is not). Note that 
n — 5fc = 1 and n — 5k = 2 are equivalent to k = (p — 5)/15 and k = {p — 8)/15, 
respectively. However, (p — 5)/15 cannot be an integer (since p ^ 5). We thus 
find that (57) vanishes modulo p unless p = 8 (mod 15) and k = [p/15j, in 
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which case (57) is congruent to —1 modulo p. Combined with the analogous 
discussion for the corresponding second term in (56), we conclude that 


{n - 5fc)3n (n - 5fc + 1)3^ 
(3n)! (3n)! 


1, if fc = [_p/15j and p = 2 (mod 15), 
— 1, if fc = [_p/15j and p = 8 (mod 15), 
0, otherwise. 


Applying this to the sum (56) and combining the signs properly, we arrive at 
the congruences (55) when p = 2 (mod 3). 

The case p = 1 (mod 3) is similar and a little bit simpler. □ 

In summary, we conjecture that the proportion of primes not dividing any 
term of the Apery-like sequences in Tables 1 and 2 is as follows. 


Conjecture 6.5. 


• Let C{n) be one of the sequences of Table 1 or sequence (77). Then the 

proportion of primes not dividing any C{n) is . 

• LetC(n) be one of the sequences (d), (a), (e), (C), (7) from Table 2. Then 
the proportion of primes not dividing any C{n) is 

In stark contrast. Cooper’s sporadic sequences sy, sio, sis from Table 2 are 
divisible by all primes. Indeed, let C{n) denote any of these three sequences. 
Then, 

C{p — 1) = 0 (mod p) 

for all primes p. In fact, we can prove much more. For any given prime p, the 
last quarter (or third) of the first p terms of these sequences are divisible by p. 
In the case of sequence sio, the sum of fourth powers of binomial coefficients, 
this is proved by N. Calkin [Cal98]. Indeed, among other divisibility results on 
sums of powers of binomials, Calkin shows that, for all integers a ^ 0, the sums 



are divisible by all primes p in the range 

n 

n<p<n + l + - -7. 

2a — 1 

In particular, in the case a = 2, we conclude that sio(n) is divisible by all primes 
p that satisfy n < p < ^ 1. Equivalently, we have 

■sio(p-j) = 0 (modp) 

whenever I ^ j ^ {p + 2)/4. Our final result proves the same phenomenon for 
Cooper’s sporadic sequences sy, sig. We note that in each case, empirically, the 
bounds on j cannot be improved (with the expection of the case p = 8 for sig; 
see Remark 6.7). 
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Theorem 6.6. For any prime p, we have 


S 7 {p — j) = 0 (mod p) 
whenever 1 ^ j ^ (p + l)/3, and 

si8ip-j) = 0 (modp) 


whenever 1 ^ j ^ (p + 2)/4. 

Proof. For the sequence S 7 , we want to show 




=0 


(mod p), 


for 1 ^ j ^ (p + l)/3. Note that for 2/c < p — j or fc > p — j the summand is 
already zero. Therefore, we assume that p — j^k^{p — j)/2. We write the 
summand as 


P - jV /P - j + fcW 2fc \ {p-j + k)\{2k)\ 

/c y V k J\p-j) kP{p-j-k)P{2k-p + jy.’ 


and observe that the denominator is not divisible by p if j 1. On the other 
hand, the factorial (p — j + fc)! in the numerator is divisible by p since 


p-j+k^p-j+ 


P-j 

2 


^ P, 


where we used J ^ (p + l)/3 to verify the final inequality. Thus, for 1 ^ j ^ 
(p + l)/3, the congruences S 7 (p — j) = 0 hold modulo p, as claimed. 

We proceed similarly for sisip — j), which is given by 


L(p-h/3J 

E (-1)* 

k=0 


V - j\ ( 2k\ (2{p - j - k) 


\k J\ p — j — k 


2 {p -j)-sk-l 
p-j 


p(p-j)-3F 

V p-j 


and, using (43), write the summand as 
(-l)'=( 2 A:)!( 2 (p-j-fc))! 


kP{p — j — k)P 


(p - j - 3fc + l)p_j_i(3p - 3j - 6k). (59) 


None of the terms in the denominator is divisible by p since j ^ 1. On the other 
hand, ( 2 (p — j — fc))! in the numerator is divisible by p since 


2 (p - j - fc) ^ 2 p - j - 


P-J 


> P, 


where we used j ^ {p + 2)/4 for the final inequality. Therefore, for 1 ^ j ^ 
(p + 2)/4, each of the terms in the sum for sig(p — j) is a multiple of p, and we 
obtain the desired congruences. □ 
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Remark 6.7. Employing (59), we observe that Si8(n) = 0 (mod 3) for n ^ 1, 
which reaffirms Corollary 5.2 for this sequence. 

Finally, as noted in [Cool2], each of the sequences in Table 1 times is 
an integer solution of (6) with d ^ 0. Observe that is divisible by a prime p 
for all n such that n < p ^ 2n. This results in a (weaker) analog of Theorem 6.6 
for these Apery-like sequences, and implies, in particular, that these sequences 
are again divisible by all prime numbers. 


7 Conclusion and open questions 

In Sections 3 and 4, we showed that all sporadic solutions of (5) and (6), given 
in Tables 1 and 2, uniformly satisfy Lucas congruences. However, for two of 
these sequences, especially sequence (r?), we had to resort to a rather technical 
analysis. We therefore wonder if there is a representation of these sequences 
from which the Lucas congruences can be deduced more naturally, based on, for 
instance the approaches of [SvS09] and [MV13], or [RY13]. More generally, it 
would be desirable to have a uniform approach to these congruences, possibly 
directly from the shape of the defining recurrences and associated differential 
equations. In another direction, it would be interesting to show that, as numer¬ 
ical evidence suggests, all of the Apery-like sequences in fact satisfy the Dwork 
congruences (3). 

The congruences (4) show that the Apery numbers are periodic modulo 8, 
alternating between the values 1 and 5. As a consequence, the other residue 
classes 0, 2, 3,4,6, 7 modulo 8 are never attained. On the other hand, the obser¬ 
vations in Section 6 show that certain primes do not divide any Apery number. 
This can be rephrased as saying that the residue class 0 is not attained by the 
Apery numbers modulo these primes. This leads us to the question of which 
residue classes are not attained by Apery-like numbers modulo prime powers p“. 
In particular, are there interesting cases which are not explained by Sections 5 
and 6? 

The second part of congruence (53) makes it explicit that every prime p = 3 
(mod 4) divides a term of the Apery-like sequence (52). Is there a similarly 
explicit result which demonstrates that the Apery numbers are divisible by 
infinitely many distinct primes? 
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